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We study electromagnetic perturbations around a Kerr black hole surrounded by a thin disk on 
the equatorial plane. Our main purpose is to reveal the black hole superradiance of electromagnetic 
waves emitted from the disk surface. The outgoing Kerr-Schild field is used to describe the disk 
emission, and the superradiant scattering is represented by a vacuum wave field which is added to 
satisfy the ingoing condition on the horizon. The formula to calculate the energy flux on the disk 
surface is presented, and the energy transport in the disk-black hole system is investigated. Within 
the low-frequency approximation we find that the energy extracted from the rotating black hole is 
mainly transported back to the disk, and the energy spectrum of electromagnetic waves observed 
at infinity is also discussed. 

o 
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I. INTRODUCTION 

It is widely believed that there exists a rotating black hole surrounded by a disk in the central region of highly 
energetic astrophysical objects, such as active galactic nuclei, x-ray binary system and gamma-ray burst sources. In 
the disk-black hole system, the propagation of electromagnetic waves emitted from the disk surface will be strongly 
affected by the gravity of the central black hole. For example, the black hole shadow is expected to be an important 
phenomenon from which the black hole spin can be observationally estimated [U-IH . The basic analysis of the shadow 
profile is relied on the ray-tracing method, and some features depending on various disk models have been investigated 
t-H ■ in 0, H[ . The further developments have been attempted in recent papers which take into account the effects of plasma 
J> ' and wave scattering [(| @] . Another key aspect of wave propagation in a rotating black hole spacetime will be the 
effect of superradiance. As was shown in [1, , the amplitude of incident waves propagating to the horizon can be 
amplified through the scattering process to extract the rotational energy of a black hole. Though this efficiency has 
been calculated for waves incident from infinitely distant regions, in the disk-black hole system it is important to 
consider the electromagnetic waves which should occur on the disk surface. In this paper we would like to assume a 
thin disk and deal with the problem of propagation of waves incident from the equatorial plane in the context of the 
superradiant scattering rather than from the viewpoint of the black hole shadow. 

The electromagnetic waves emitted from the disk surface should be directly carried away to infinity, and partly 
absorbed by the black hole. Our main interest is focused on the electromagnetic energy transport in the disk-black hole 
system. We will clarify how the wave absorption across the horizon generates superradiant energy outflows from the 
black hole to disk and infinity. This energy transport to the disk may be interpreted as a feedback mechanism which 
plays a role of disk reheating, while the contribution to the energy flux at infinity may be useful for an observational 
check of the superradiance effect. (Such energy outflow may also contributes to a change of the black hole shadow. 
However, to pursue this possibility is beyond the scope of this paper.) 

Because the ray-tracing method is not applicable to the analysis of the superradiant scattering of electromagnetic 
waves, our approach to the problem starts from solving the vacuum Maxwell equations in Kerr geometry. The assumed 
boundary condition is the existence of a thin disk (i.e., a surface current) on the equatorial plane. Namely, it is required 
that some components of electromagnetic fields become discontinuous at the equatorial plane, and outgoing energy 
fluxes are emitted from both the upper and lower sides of the disk surface. 

The Kerr-Schild (K-S) formalism for solving the Einstein-Maxwell equations is a useful method to overcome the 
mathematical difficulty due to the disk boundary contribution |lfj,|ll|. In fact, if the electromagnetic fields are treated 
as p erturbations in Kerr geometry, all the field components are simply derived by two arbitrary complex functions 
12], which can be appropriately chosen according to the disk boundary condition. Hence, in SecHH we introduce the 
outgoing K-S field as a model of the disk emission and discuss a singular behavior of the two complex functions 
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at the equatorial plane. Unfortunately, the outgoing K-S field fails to satisfy the condition that no outgoing waves 
are present on the horizon. Hence, the physical field satisfying the horizon boundary condition should be modified to 
the form Fp V — Fj^f + F™ , where the additional vacuum field be interpreted as the non-Kerr part due to the 
supcrradiant scattering of waves emitted from the disk and can be continuous even at the equatorial plane. In order 
to facilitate the analysis of the scattering problem, we consider the Newman- Penrose quantities 0, HH corresponding 
to the electromagnetic field F M „ = F^ + F™ , and in Sec. IIIII we express them as the infinite sums of a complete 
set of modes. Based on the mode decomposition, in Sec. IIV1 we derive the formulae to calculate the energy fluxes 
at the boundary surfaces including the horizon, the equatorial disk and infinity. In Sec. |V1 the superradiant energy 
transport from the black hole to the disk and infinity is explicitly estimated within the low-frequency limit of the 
wave fields. Hereafter we use units such that c = G = 1. 



II. KERR-SCHILD FIELD AND NEWMAN-PENROSE QUANTITIES 

Let us consider the electromagnetic waves emitted from disk surface around Kerr black hole, using the framework of 
the Kerr-Schild formalism (see details in [10]). Though this formalism is introduced to solve the full Einstein-Maxwell 
equations, it may be applied to obtain electromagnetic perturbations on Kerr background. The metrical ansatz is 

gil u = _ 2He 3fi e au , (1) 

where , rf iv is the metric of an auxiliary Minkowski spacetime, H is scalar function, and e 3AI is a null vector field, which 
is tangent to a geodesic and shear-free principal null congruence. It is convenient to calculate the Einstein-Maxwell 
equations using tetrad components. All other null tetrad vectors are defined by the condition 
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= 3 ab , (2) 



where latin and greek suffixes mean tetrad and tensor suffixes, respectively. A tensor T""' is related to its tetrad 
components T a b y by either of the two equivalent relations 

= <:;<!: - t* Z) t/ ::: = ... r„' z . (3) 

The essential point of the Kerr-Schild formalism is to use the complex form of electromagnetic field tensors given 

by 



F av — F^-y ~\- iCfivpcrF^ , (4) 



where e^pa is completely skew-symmetric, and equal to £1234 = (— g) 1 ^ 2 - The corresponding null tetrad components 



are 

^ab = Fab + -ji^abcdF Cd , (5) 

where e a b c d is completely skew-symmetric, and £1234 = i. By virtue of the definition (JS|) and the Einstein equations, 
the tetrad components F32, F41, and F42 are found to be zero. The electromagnetic fields are completely determined 
by only two complex components F12, and F31. It is interesting to note that the Kerr-Schild form remains valid, even 
if a back reaction on the gravitational field by the electromagnetic field is considered. 
A part of the Maxwell equations allows to write the tetrad components as 

^"12-^34 = AZ 2 , (6) 
J- 3 i = 1 Z-{AZ), U (7) 

where Z is the complex expansion of the null vector e 3 , and commas denote the directional derivatives along chosen 
null tetrad vectors. The functions A and 7 should be determined by solving the other Maxwell equations. 

In this paper we treat the electromagnetic fields as perturbations on Kerr background, and use the outgoing Kerr- 
Schild coordinate system to describe electromagnetic waves emitted from disk to infinity. The outgoing Kerr-Schild 
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form of the Kerr metric is given by 

ds 2 = -dP +dr 2 + Sd6» 2 + (r 2 +a 2 )sin 2 ^" 2 



+2a sin 2 0drd<^ 
2 Mr 

H — — (di-dr -asin 2 0d(p) 2 , (8) 

where £ = r 2 + a 2 cos 2 0, and M and a denote the mass and the angular momentum per unit mass of the black hole, 
respectively. The function H in Eq. ([!} is given by 2Mr/? 2 S, where P — l/\/2cos 2 (0/2), and the null tetrad vectors 
are given by 

e 1 = [-tan(0/2), tan(0/2), (r - ia cos 0), (a + ir) sin 0] , (9a) 

e 2 = 2~h-^ [-tan(0/2), tan(0/2), (r + iacos9), (a-ir)sin0], (9b) 



e 



P [-1, 1, 0, asin 2 0] , (9c) 
e 4 = 2-*[l, cos0, -rsin0, 0] + He 3 , (9d) 

with e 3 an outgoing null geodesic. Then, the functions A and 7 for electromagnetic perturbations on Kerr background, 
can be written as 

A = *™. (10, 

where r = t — r + iacosO because e 3 is chosen as an outgoing vector field, Y = e 1 ^ tan(0/2) (see [13): and comma 
means differentiation with respect to a given variable. Further, we obtain Z/P — l/(r — iacos0) for the complex 
expansion Z in Eqs. (J6)) and ([TJ. It should be noted that the tetrad components T a b can written by the two arbitrary 
complex functions ip(Y, r) and 4>(Y.t) as follows, 



(r — iacos( 



= ,...„ / ,, 2 , (12) 



ia cos \ I 



. . ~2cos 2 (0/2) r - ia , . n r + ia 
— 1 tan(0/2) r-,acos0 ^ +Smg (,-, a cos0) 2 ^ 



+ 0(F,r) ^-^V (13) 

r — ia cos J 

Hereafter we call this solution the Kerr-Schild field. 

To see clearly superradiant energy transport in the disk-black hole system, it is convenient to introduce the Boyer- 
Lindquist coordinates, which lead to the metric 

, / 2Mr\ „ AaMr , £ , , A , , , , 

ds 2 = -II dt 2 sin 2 66.t6.ip + —dr 2 + Sd0 2 + - sin 2 0d<^ 2 , (14) 

V S / s s 

where A = r 2 + a 2 — 2Mr, and A = (r 2 + a 2 ) 2 — a 2 A sin 2 0. The Boyer-Lindquist coordinates t and ip are related to 
the outgoing Kerr-Schild coordinates t and (p as follows, 

- 2Mr , , a , 

df = dt + — ^— dr, dy3 = dip + —dr. (15) 

If the electromagnetic perturbations written in the Boyer-Lindquist coordinate system are assumed to be 
functions of three variables at — ip, r, only, we can simply expect that the superradiant scatte ring s occurs under the 
condition a < fin for the angular velocity J7h of the black hole and a frequency parameter a [111 ]. The expectation 
motivates us to specify the complex function tp(Y,r) and 4>(Y, r) in Eqs. (TT2"|) and (|13l) to the forms 

ip(Y, t) = iJ,(X), 4>{Y, t) = {-iaaX^,^ x + *{X))/Y. (16) 

where the complex variable X is defined by 

X = e- taT Y = e - iCTT +^ tan(0/2), (17) 
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and the term —iaaXip.x is included in <p to simplify the expression of Ff,,, which will be given later. It is easy to see 
from Eqs. (TT2l and (fTBf that by virtue of the choice of ip and <p the field components F\ 2 and J-31 depend on t and ip 
via the variable err — (p in X, and we obtain 



<jt — (p = at — if + iaa cos 9 — or* + fl^iL(r) 
where the tortoise coordinate r* is defined as 



r* = 



A 



-dr = r — n 



2Mr x 
ri - r 2 



In 



' 1 



n - r 2 



2Mr 2 
r 2 - ri 



In 



r-2 



n - r 2 



using the outer and inner horizon radii ri and r 2 , respectively. For the function L(r) given by 



L(r) 



2Mn 

A 



dr 



2Mr x 
r\ - r 2 



In 



n 



r - r 2 



(18) 



(19) 



(20) 



we can check the asymptotic behaviors such that L — r* — > on the outer horizon r = ri, and L — > at infinity 
r — > 00. 

It is a straight forward task to derive the field components in the Boyer-Lindquist coordinate system from J-12 and 
F31 given by Eqs. (fT2|) and (jlU)) . Using the specified form of ^ and 0, the K-S field components denoted by are 
given by 



f ks 

r tr 



r t<p — 



-Re 
-Re 
-Re 
Re 



(r — ia cos 9) 2 A sin 

. . ~ m 2 



ia sin ( 



(r — ietcos( 



sm 



(6 + zS) 

(e + zs) 



sin 9 

.2 1 „2 



(9 + zS) 



— z(r + a ) sin 



KS 

rip 

KS 



= Re 



—a sm 



(r — iacosfl) 2 
4>{X) ,r 2 H 



2E 
Asin6»" 



-Re 



(r — iacos^) 2 
1 



A sin 6» 



+ 2a (9 + i3) 
(9 + zS) 



(9 + zS) 



where the function 9 and S are defined as 



9 = $>(X) + Xifj iX {2a<r cos 2 (9/ 2) - l) /(r - ia cos 9), 
S = X ^ yX acr[2r cos 2 {9/ 2) / (r - ia cos 9) - l] . 



(21a) 
(21b) 
(21c) 
(21d) 
(21e) 
(21f) 



(22) 
(23) 



Here we consider the boundary condition on the disk located at the equatorial plane 9 = n/2. It is well-known that 
any complex function which is not a constant should have a singularity on the complex plane. We will assume that the 
existence of a singularity in ip(X) and *&{X) on the complex X-plane, is due to a surface current on the equatorial plane 
9 = 7r/2. This means that the components F^ s , F^, and F^° (namely, the imaginary part of ip and the real part of 9) 
become discontinuous at 9 = n/2. Such a discontinuity will be generated if a branch point in -0 exists at X = e 1 ^ where 
f3 is a real constant. For example, as was discussed in [ll|, the function ip(X) = ipo [(X~ 2 + X 2 ) 3 / 2 — X~ 3 — X 3 ~\ 2 
with a real constant ipo has four branch points at X 2 = ±i, and the imaginary part of ip and the real part of Xip^x 
become discontinuous at 9 = ir/2. In this case the ratio 1 $>(X)/Xip,x may be chosen to be a real constant, as will be 



done in (1691. Further it should be noted that the absolute value \X\ 



~aa cos 9 



9 = n/2. The branch point X 



tan(0/2) becomes equal to unity at 



may also appear on some conical plane 9 = 9q(^ n/2) giving \X\ = 1 if aa > 1. 



Hence, in the following, the allowed range of the frequency parameter a is limited to the range < a < I /a, for which 
we obtain \X\ < 1 in the upper region < 9 < n/2 and \X\ > 1 in the lower region 7r/2 < 9 < n. 

We must also consider the regularity condition for F^f at the polar axis (i.e., at 9 — 0, n). Noting that \X\ ~ 
sin(0/2) in the limit 9 — > and \X\ ~ l/cos(0/2) in the limit 9 — > n, we find the boundary condition for ip and \& to 
be for 6^0, 



iP(X)~^{X)~X 2 , 



(24) 
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and for — > n, 

ip(X) ~ *(X) - l/X 2 . (25) 

Finally, let us discuss the boundary condition on the horizon and at infinity, by introducing the electromagnetic 
Newman-Penrose quantities (b a (a — 0, 1, 2) [l3j]. Using the Kinnersley's tetrad well-behaved on the past horizon 
such that 

i" = [(r 2 + a 2 )/A, 1, 0, a/A], (26a) 
n^ 1 = [r 2 + a 2 , -A, 0, o]/2E, (26b) 



m M = [ia sin 



0, 0, 1, i/sin0]/2 1/2 (r + iacos0), (26c) 



the Newman-Penrose quantities are defined as 



b = F^m u , (27a) 
1 

r 

b 2 = F^m»n v , (27c) 



-F,„(iV+mV), (27b) 



and satisfy the Maxwell equations written by 

(r — iacos6) ( iasin 



V2 

1 



V2(r — ia cos 
A 



\/2{r — ia cos 



/ „ ia sin \ r , 

£ i~7 = « Uo = Z>o [tfi (r-ia cos 0) 2 ] 28a 

\ (r — ia cos 0) / L J 

£ [0i (r - ia cos 0) 2 ] = X> [02 (r - ia cos 0)] , (28b) 
X>J [0i (r - ia cos 0) 2 ] = £| [0 2 (r - ia cos 0)] , (28c) 



(r — iacos0)A /_+ 1 \ , „i r , , „ l01 

1=—*- V\ - ; Uo = 4 [Mr - *«cos0 2 , 28d 

V2 V ( r — «acos0) / 

where the differential operators T> n and £„ are defined by 

V n = l^dp, + 2n(r — M)/A, (29a) 

2?t = -2(£/A)n' i 9 A( +2n(r-M)/A, (29b) 

£„ = \/2(r — ia cos 9) fh^d^ + ncot6, (29c) 

£t = \/2(r + ia cos 0)m^ + n cot 0. (29d) 
From Eqs. (|2Tj) . (|26|) and (|27|) we obtain the Newman-Penrose quantities for the Kerr-Schild field as follows, 

$ S = 0, (30a) 

- ^ X L™ > (30b) 



2(r — ia cos 0) 2 ' 

r 5fi/i „ I 

(30c) 



^ S = -p " 

V2(r — ia cos 0) sin 



+ / ^ ( icrr cos + ao — 1) 

(r — ia cos 0) 



Because the electromagnetic waves are assumed to be emitted from the disk, no ingoing waves should exist at infinity. 
Hence, we require the asymptotic behavior of the Newman-Penrose quantities obeying Eqs. (|28[) to be 

0o~l/r 3 , 2 ~(l/r) 5 (i-r,), (31) 

in the limit r — > oo. On the other hand no outgoing waves should not exist on the horizon, and we require 

0o — (l/A)/(i + r„), 02 ~ A, (32) 

in the limit A — > 0. It is easy to see that the Kerr-Schild field satisfies the boundary condition only at infinity. The 
horizon boundary condition breaks down, because 0?P does not vanish at A = 0. Therefore, to obtain the physical 
field a which is well-behaved on the horizon, some vacuum field denoted by 0„ is added to the Kerr-Schild field 
0?- s as follows 



aSW , _ ,KS , iSW , _ iKS , jSW 

»o ' 01 — 01 + 9i j 92 — <Pi +92 > 



(33) 

where 02 is required to vanish on the horizon. Though the Kerr-Schild field describes the disk emission, the additional 
field is expected to represent the effect of wave scattering (or absorption) by the black hole. In the next section we 
will describe the scheme to obtain the additional field 0^ W , by imposing the conditions (|3ip and (|3"2"j) on a . 
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III. WAVE SCATTERING 



As the first step to analyze the scattered wave field (f>a W , let us expand the functions ip(X) and ^(X) in Eqs. (|30|) 

as 

i>(x) = Y, a ™ x " l > ( 34 ) 

m 

*{X) = Y,b m X m , (35) 

m 

where from the condition (|24[) m runs from 2 to oo for < \X\ < 1 (corresponding to the upper region < 8 < ir/2), 
while from the condition ([25]) it runs from —2 to — oo for 1 < \X\ < oo (corresponding to the lower region ir/2 < 9 < it). 
Such an expansion will be possible, because tj}(X) and ^{X) are assumed to be regular except at branch points on 
the equatorial plane \X\ = 1. Note that the m-th terms ip m and ^ m are proportional to exp [im(at — ip)], which 
represents a mode with the wave frequency 

oj m = ma, (36) 

for m > 0, while the wave frequency should be understood to be — u m for m < 0. By virtue of the expansion of ip 
and the Kerr-Schild field </>!p is rewritten in to the form 

oo 

(r- wcos6l) 2 0!p = e- luJmt+lm,p S^ s (r,e)e i{uJmr '' mnilL) , (37) 

m— — OO 

pCbOJm cos j.„„m 

S* s (r,9) = H m = [(r - iacos6)b m + a m (iuj m r cos 9 + auj m - m)] . 

V 2 sin 

(38) 

Because the modes for m = ±1, should not be included in Eq. (|37l) . we have H±\ = Hq = 0. Further the expansion 
forms (I34[) and (1351) mean that for m > 2 (or m < —2) ij) and $ must vanish in the lower (or upper) region. Hence, 
the factor H m in (|38[) is given by the step function such that H m = (1 + m/\m\)/2 in the range < # < 7r/2 and 
ii,,, = (1 — m/|m|)/2 in the range tt/2 < 9 < tt. From Eq. (|37|) we have the asymptotic behavior near the horizon 
r = 7"i as follows, 

(r-^acos0) 2 < ^ ! ) P-Xl e ~ ^ " m^+^mV ' 5 '™ S ( ^1 ' 6 ') e^ ' ^mr *' fc ™ = w ™ _ m ° H ( 39 ) 

m 

which should be canceled out by the scattered- wave field ^>§ W according to the horizon boundary condition. The easier 
way to construct such a vacuum non-Kerr-Schild field will be to use the expansion form written by the spin- weighted 
angular functions Sj m The application of this mode decomposition to 4>q W and (/>f w leads to the result 

su _ ^e-^+^S^mttir), (40) 



u 



m.l 



{r- l acos9Y4>l w = ^e^^S^^^ir), (41) 

„l 



where R$ is the radial function and I > \m\ > 2. Note that the component <j)f w can be derived by using the Maxwell 
equations (|2"51) . Therefore we consider hereafter only the two components </>fj w and </>f W - 

The asymptotic behaviors of R^ near the horizon and at infinity are well-known. For example, in the limit r — > n, 
we give the radial function as follows 

i?« ~ C^A-'e-^+C^e^, (42) 
R[ m l) ~ D^Ae-^+^e^, (43) 

where the outgoing parts with the amplitudes C^* and D^ t are also included to satisfy the condition </>f w + 4>f s — > 
on the horizon. On the other hand, in the limit r — > oo where no incoming waves exist we obtain 

* E^e-^'/r 3 , (44) 
Rt 1] * F^re-^*. (45) 
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The coefficient ratios Cj^/D 1 ^ and E^/F^ have been derived in [9] using the Teukolsky equations, and the results 
are given by 

r m /n in 32ik m M 2 r 2 (-ik m + 2e) 

E?^/F^ = -JL (47) 

ZLO m 

where B = (E + a 2 cj„ 2 — 2auj m m) 2 + 4maoj m — 4a 2 w„ 2 and e = (ri — ra)/ AMr\ , and -E is the eigenvalue of the angular 
equation. In the low-frequency limit aw m — » 0, we have E — > 1(1 + 1), which is the case analyzed in Sec. [V] Further, 
we obtain the ratio Cffi/Dffi written as 

Qm /E>lm — TTn 27^1 , o \ ' V^8) 

8ik m M z r 1 ^{ik m + 2e) 

From Eqs. (|39|) and (|43|) - the asymptotic behavior of 02 near the horizon is written as 

(r - l acos0)</> 2 ~ ^ e -*-m*+fm V [s}- 1 )^) Ae"*-*-* + D°f e^*) + S^ s ( ri ,9)e ik ^] (49) 

Hence, the horizon boundary condition for 02 leads to the relation 

E ^ ut C X) W = S^iruO), (50) 

i>|m| 

from which the coefficient -Dj^* is determined by 

Am* = - S^ s (ruO)S i lm 1) (9) siaede. (51) 
Jo 

for the given Kerr-Schild field. 

The important problem to be solved in relation to the superradiant scattering of disk emission is to estimate the 
ratios Cf^/D^ and Fg^/Dft?, based on the radial equation 



A d 2 i?( 
A 



dr 2 



K 2 + 2i(r — M)K 

Airuj m — X 



= 0, (52) 



where K = (r 2 + a 2 )u> m — am, A is separation constant written by A = E — 2amus m + a^uj^. 

Finally, we summarize our proposed approach which is the derivation of the scattered wave Fjj^. In our approach, 
the property of the disk emission is given by the Kerr-Schild field (jJ£ s , that is, the two any complex functions ip and 
^ or the expansion coefficients a m and b m in Eqs. and (j3"5")) . However the horizon boundary condition (|32p breaks 
down, because ^>JP dose not vanish at A = 0. Therefore, to obtain the physical field 4> a which is well-behaved on the 
horizon, we introduced the vacuum field </>„ w which represents the effect of wave scattering (or absorption) by the 
black hole. Then, the coefficient D^ t in the outgoing part of the scattered wave should be determined by Eq. (|5"T|) . 
using the Kerr-Schild field <^ s . To determine the scattered field 0q W , we must solve the radial equation (j52")l . using 
the boundary value given by Eq. (IBTl) on the horizon. If the radial equation ([52]) is solved, the ratios D^/D^ 
and F^^/D^ 1 will be obtained. Therefore, the scattered wave is represented only by the coefficients D^ t 
connected with Kerr-Schild field on the horizon. 

Before pursuing the analysis of Eq. (|52l) in details, we must present the formulae to calculate the energy fluxes 
from the disk, on the horizon and at infinity, because our main purpose is to clarify the energy transport via the 
superradiant scattering process in the disk-black hole system. This will be done in the next section. 



IV. ENERGY FLUX 



Using the Newman- Penrose quantities obtained in the previous section, let us present the useful expressions of the 
energy flux vector defined by 

= -Tl, (53) 
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where 

Tpv = [<t>o(t)Qn^n u + 2 02^L + 20i0i (l(i n j) + 

-40 o 0in(imj) - 40102^77^) + 20200^^] + C.C.. (54) 

Because the component 0i given by 0o and 02 through the Maxwell equations (|28[) . the energy flux vector £ M can be 
written by 0o and 02 only. 

Note that the energy flux vector for the wave fields considered here is oscillatory with respect to the time t (as 
well as the azimuthal angle tp). To estimate the efficiency of the energy transport, we must consider the time-average 
quantities such that 



2tt Jo 



Adt, (55) 



with the frequency parameter a. Because it is easy to see that the time-dependence on the energy flux vectors £ M 
arises from the electromagnetic field components Q , we consider the time-averaged quantities (0 a o ) written as 

a r 2 * /a - 

(0a06> = r- / a 0fcdr;. (56) 
2tt Jo 

Note that the electromagnetic field components a are expanded as 



ba = E 0am, (57) 



m— — oo 



with 1 771 1 > 2. In particular, 0o m and 02 m are written as 



e 



-iiUrnt+imtp 



E S^m^ir), (58) 



l>\m\ 



(r — ia cos 9) 2 (f)2m 



—iujmt+imtp 



i>\ 



(59) 



From Eqs. (|58p and ([59^1 , it is easy to see that the time-averaged quantities (0 a 0b) are given by the mode decompo- 
sitionas follows, 

oo 

(0a0b) = E ^m0hm, (60) 

771 — — OO 

with | m | > 2, which will lead to the time-averaged energy flux vectors written as 

oo 

<^}= E Wn), (61) 
m— — oo 

with | m | > 2. Hereafter, we consider the mode-decomposed and time-averaged energy flux vector 

We consider the angular component (£ ) of the energy flux vector as the emission from the disk surface. Nothing 
that \/T, = r at 8 — tt/2, the energy flux £y> per unit area can be evaluated as 

£»(r)=r(S+)+r(£-), (62) 

where (£^) are equal to ±(£^) in the limit 9 — > 7r/2± corresponding to the disk emission from the upper and lower 
side, respectively. Then, are obtained as 

<£ m ) = -5-27-^ \ d ^-TT^ \ iB ™> ( 63 ) 

8irr z (m — aui m ) 4ir(m — au) m ) 

4n = [^L^rn+^L^+^^f^ + ^^^Y (64) 
B t = (00m 02m - 4>OrJ>2m)- (65) 
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To evaluate the total flux radiated from disk surface, it is easy to obtain the total energy flux E® as follows 

/>oo 

E» = 2tt / 8%rdr. (66) 



After a tedious calculation, we have the total energy flux E^ given by 

E m = T7 ~ r [A+ - O" - /'°° - B")r 2 dr. (67) 

4(m-aw m ) L Jri 2(m- aui m ) J ri 

Noting that for the Kerr-Schild field O <S = and 0Jp — > in the limit — > ir/2 — 0, the corresponding total energy 
flux [E^] KS is obtained as 

[*£] KS = 9(m 1 — > [(- 2 ^f^ s ) H - (^r<s s u . (68) 

To assume the continuity of 0^ w at = 7r/2 on the horizon for to > 2, we impose the same continuity condition cf>^ s 
which is interpreted as the condition that the disk does not extend to the horizon, and choose (?* 2 02"m S< ^2m )h to be 
zero. Then, from Eq. (|39|) it is easy to see that the coefficient b m is determined as follows, 

to - auj m 

b m = a m . (69) 

n 

In particular, noting the continuity of the scattered field at 9 — n/2, namely, 0q SW = 4>~ sw , the total energy flux 
i?Dm is obtained as follows, 



E,, 



2(to — aUJm) 
(r 2 6 +KS d> +sw ) 

V < V2m < V2m )°o 

2(to — auj m ) 



Im(«^ s ^ w )r 2 dr. (70) 



For m < — 2, we can use the same formula only by exchanging the subscript + for — and vice versa. 

Next let us calculate the energy flux on the horizon, where we obtain the radial component (£ m ) of the energy flux 
vector as 



87rfc m SH 



... „ / (A 2 Om 0Om)H . 2 , 
-(A 00m 9>0m JH ~ r, . Q <J0 ^ Sm I 



(71) 



_2MV TU '"" """" 2sm9"\ S H 
Further, we can evaluate the total flux E^ integrated over the whole horizon surface as follows, 

E* = 2tt / (££) H E H siii0d0, (72) 
Jo 

which reduces to the form 

E m = - 8fc ^n [ ( A2 ^m0om) H singdg. (73) 

Note that the integration of the second term in Eq. ([7T|) is canceled out, because 0o m is continuous even at 9 — ir/2. 
The result given by Eq. (I73[) shows that the energy extraction from the black hole occurs for incident waves with the 
frequency parameter a in the range < \u) m \ < |to|Oh (i.e., k m < for to < and k m > for to > 0) which means 
the range < a < fin, in accordance with the result of the usual superradiant scattering [§]. We can rewrite the net 
flux £" m into the form 

8k m Mri ^ 



i 



in which the coefficient C/ m will be given by solving the radial equation ([52]) under the low- frequency limit (see Sec.|V|. 

Finally we turn our attention to the radial component of energy flux vector (£ r ) m at infinity. The energy flux 
vector is written by 

(£ T m )oo = ;T^(r 2 02m02m)oo- (75) 

Znr 
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It is easy to see that no ingoing energy flux exists at infinity. Further we calculate the total flux at infinity as follows, 

E% = 2ir f (C)oor 2 sin 6d6. (76) 
Jo 

Then the total flux at infinity is given by 

= I (r 2 <f> 2m fc m )sm9d9. (77) 



Considering the contributions form the Kerr-Schild field and the scattered wave, the total flux at infinity is rewritten 
by 

B%= rr^\& + 2Re(&^) + \4^] oo sined9, (78) 
Jo 

where from Eqs. (|30c|) , (|35|) . and (p9")) we obtain the asymptotic form of Kerr-Schild field as 



C -e-**" *H m ame - tm n WZ) [(m - au m ) - iu> m r x costf] , 

V2rri sm t) 



(79) 

V z,rr\ sm a 

while nas the asymptotic form given by Eq. (|4"5j). 

Here, we assume the Kerr-Schild field to be reflection-symmetric with respect to the equatorial plane. This symmetry 
is corresponding to the condition that the two expansion coefficients are written by 

Q, m Q— rni h rn 6_ m , (^0) 

which allows us to calculate the energy flux for m < 0, by using the result for m > 0. 

In following section, we see the energy transport to calculate the total energy in each region using the low-frequency 
limit. Then we should solve the radial equation (I52[) to determine the amplitude of the electromagnetic fields on the 
horizon and at infinity. 



V. LOW-FREQUENCY LIMIT 



In the previous section, we discussed the time-averaged energy flux vector on the boundary surface. To evaluate 
explicitly the efficiency of the superradiant scattering, we attempt to solve the radial equation (|52|) using the low- 
frequency limit (au! m — > 0) according to the procedure in Q. 

For auj m <C 1, the angular function Sj^(9) is known to be given by 




21 + 1 (l + m)\ {I - m)\ . 2l 



I- s 
r 



{1 + 8)1 (l-s)l 

I + s 

r + s — m 



sin" {9/2) 



{-l) l -' r - s cot 2r+s - m {9/2) 



(81) 



where the eigenvalue E is obtained by E — {I — s){l + s + 1) (see [1J|)- On the other hand, we obtain the radial 
function through the asymptotic matching method. To consider the dimenssionless radial equation from Eq. (|52[) , we 
introduce the dimenssionless condition and parameter as 



r — ri 
n — ti 



Q 



n - r 2 

Then, for u) m M <C 1 and x S> max(Q, I) (i.e., in the distant region ), the radial equation can be rewritten by 

., _> 2iu m {ri - r 2 ) 1(1 + 1)~ 



d 2 i?!- 1} 



Ira 



dx 2 



w m (ri - r 2 y 



- 



(82) 
(83) 

(84) 
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This equation is expressible in terms of the confluent hypergeometric functions, and we obtain the outgoing wave 
solution given by 



2nxe~ KX r°° f , , / t 



1+1 



with k = —i{ri — r-i)ijj mi and C = F^ t /(2iui m ). It is easy to see that the asymptotic behavior of Eq. (|85|) is written 
by written by 

R\ m 1] ~ C2Kxe-"*, (86) 



for \kx\ ^> 1 in accordance with Eq. (|45|) . and 



? (-D „ n ( 2 ! 

(/-!)!' 



i?^^C-^-(2 K )-'- 1 x-', (87) 



for <C 1. In the region near the horizon (x <C l/ui m (ri — r^)) the radial equation (|52j) can be rewritten as 

d 2n(-l) 

[x(x + l)] 2 ^ + [Q 2 - »Q(1 + 2i) - + l)z(x + 1)] fl, ( ~ 13 = 0, 

where it is obtained by neglecting in Eq. (1521) all the terms containing uj m except the one which enters into Q, and 
the solution with ingoing and outgoing boundary condition at the horizon is given by 



RL^ = + 4^ (Out), (89) 



where i?[ m l3 (In) and ii^ 13 (Out) are written by 



aW M - (2/ + i)! r(-MQ) r ( i + 2t Q) / g V Q 

Klm{ln) ~ A (l- 1)! r(2 + 2iQ) r(I + 1 - 2*Q) U + 1 J + lj 



T(2 + 2 ? Q) U r(n+l-f)r(n + f-l)(-o;)" 
r(l-na+l)!^ r(n + 2 + 2iQ) n! ' 1 ' 



^ m (Out) - + r(/ + 1 + 2iQ) [ x+1 



T(2tQ) ' r(n-l-/)r(n + Z)(-a;)" 

r(-l-I)(i-l)!^J I> + 2£Q) n! ' 1 j 



respectively. From Eqs. ([50]) and ([?!]) , the amplitudes Dfe and are given by 



lm Im 

.„ _ . ,„(2i + l)! r(-2»Q) T(l + 2zQ) 



A™ - -(n-r- 2 ) ^ (z_i)!'r(2 + 2zg)r(z + i-2 4 Q)' (92) 

put _ ^ (2/ + !)! r(l + 2zQ) 

/m ~ (l + l)! r(I + l + 2»Q)' 1 ' 

We obtain the solution with the asymptotic behavior as follows, 

i?! m 1} ~ {-l)- l Ax-\ (94) 

for x S> 1. If all the parameters satisfy the condition max(Q, I) <C ljio m (r\ — r%), we have the overlap region in which 
both the outer expression (|8T[) and the inner expression (|94|) hold. In this region, we can match the leading-terms of 
the solutions (|571) and (fMj) . and this matching yields 

^ = (-l) ; (2fc)- i - 1 ^^ y C. (95) 
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Then we can obtain the ratios F^/DP"* and Dj^/D^ as follows, 



pOut 
r Ira 
nOut 

nOut 



(-l)- l (r 1 ~r 2 )-\2 K ) 



l+1 (l-l)\(l + l)\T(l + l + 2iQ) 



(2/)!(2/ + l)! T(l + 2iQ) 
(/ + !)! r{l + l + 2iQ)T(-2iQ) 



(n - r 2 ) 2 (I - l)! r(Z + 1 - 2?:Q)r(2 + 2iQ) 



(96) 
(97) 



Here the coefficients -F)„ ut and Dj^ represent the amplitude of the outgoing waves at infinity and the ingoing waves 
on the horizon in the component </>2 W , respectively. Further, from Eqs. and (|97l) it is easy to see the coefficient 
C& describing the ingoing waves in the component 0q W to be 



°2m 



Note that from Eq. the ratio \C\ 



In 12 



D 



TOT 



Out 



,r(/ + l + 2zQ) 

J r(; + i-2iQ)' 



is given by 



\C, 



}m\ 2 / 



\D 



Out 1 2 
Im I 



= 4. 



(98) 



(99) 



Next we pay attention to the energy extraction from the black hole. From Eq. (|7Tj) . the total flux on the horizon is 
given by the component <f> . On the other hand, the scattered field is induced and connected by the Kerr-Schild field 
on the horizon through Eq. (|51|) . Further, using Eq. (|69j) which requires the absence of the disk on the horizon, the 



coefficient Dj^ 1 ' is obtained as 



nOut _ 



Q,k m 

ns2 H 



Therefore, using the ratio Cfa/Dft? given by Eq. (|98p we can evaluate the component cf>o near the horizon as 

,(i) , , M ' -f- 1 



»0m 



F 



— imuj 1Tl t-\-irmp q{±) a - 



2iQ) nOut 

r(z + 1 - 2«q) /m ' 



(100) 



(101) 



which is useful to see the distribution of the energy flux on the horizon through Eq. (fTTj) . 

The ^-dependence of A 2 |0om| 2 on the horizon is shown in Fig. [I] from which we can see that a highly asymmetric 



[A 2 |0, 




FIG. 1: The amplitude A 2 \4>0m\ 2 of the mode with m = 2 on the horizon as a function of zenithal angle 8. The thin line, 
dashed line and heavy line correspond to the cases Q = (a = 0), Q = 10 (a ~ 0.995), and Q — > oo (a — > 1), respectively. 



profile appears as the spin parameter a increases. For example, for the modes with m > 2 the peak of A 2 |0om| 2 
exists in the lower hemisphere ir/2 < 8 < tt, though the K-S field (j^m vanishes in this range. On the contrary, 
for the modes with m < —2, the peak of A 2 |0om| 2 ° n the horizon exists in the upper hemisphere < 8 < ir/2 
according to the relation A 2 \4>o m (6)\ 2 = A 2 \<j>o_ m (ir — 8/2) | 2 . Even if the net energy flux estimated by the sum 
(£ r (8))n = ^2 m>2 (£m(Q))ii + Em<-2^m(^))H nas a ^-dependence symmetric with respect to the equational plane, 
the asymmetric profile of J2 m >2(^m(^))ii given by A 2 |0 O m(#)| 2 is an interesting feature of the Kerr black hole. 

Now let us evaluate in Eq. (|BrJ)) and in Eq. (1751 under the low-frequency approximation, keeping the terms 
up to the first order in Muj m . The total flux radiated form the disk surface is given as 

D (m-ow™) 2 

h rn — 1 \ a m\ , (102) 

4n 
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while the total flux at infinity is obtained by 



The difference 



El 



E,, 



(to — 2auj m )\a r , 
4^ 



E° 



auj m 

2^7 



(103) 



(104) 



means the energy inflow from the disk to the black hole. On the other hand, from Eqs. (|74l . (|101[) and (IM|) the net 
extracted energy on the horizon is calculated as 



El 



aui r , 



(105) 



Now we can discuss the energy transported in disk-black hole system under the low-frequenct approximation. From 
Eq. (|104l) a part of disk emission turns out to be transported to the black hole. On the other hand, from Eq. (I105p . 
it is easy to see that the net flux on the horizon indicates the energy extraction from the black hole. Considering the 
energy conservation in the disk-black hole system, the energy flow (|104[) transported from disk to black hole induces 
the black hole superradiance (|105|) . and the total flux E^ + (E^ — E™) returns to the disk surface (see Fig. [5]). Then, 




FIG. 2: Schematic diagram of energy transport in the disk-black hole system. The energy radiated from disk surface is shown 
as solid arrows. The superradiant energy flows from the black hole to the disk surface and to infinity are shown as heavy dashed 
and thin dashed arrows, respectively. 



it is interesting to note that the feedback energy flow is estimated to be 



El + (EV-E%)=*-{E»-E%), (106) 
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which is one and half times as large as the energy flow — from disk to black hole. On the other hand, if we 
consider the limit a — > (i.e. the case of the Schwarzschild black hole), the net flux on the horizon is given by 

E*~-^\a m \ 2 , (107) 

that is, no energy feedback from the black hole to the disk occurs. 

Finally let us discuss the energy flux observed at infinity. From Eq. (|78|) . the net flux at infinity can be divided as 
follows, 



E m= F m+ F m + F m> ( 108 ) 



where 



„2aoj m cos 8 2m in /q\ 

*E = I H^\am\ 2 - , ™ [(m - aa, m ) 2 + a, m V cos 2 g] dfl, (109) 

Z/V ~y Sill (7 




^aw^ cos t 



= 2Re|jf H m ame ~ ^ nm(d/2) [(m - au, m ) - iu m n cos 8} £^^(6)^ , 

(110) 

F m = E [ V F^F^S^ 1 \8)S^ ) (9)sm9d6. (Ill) 

Here from Eqs. (IM|) and (|100[) the coefficient i 7 ^ 11 * is obtained as 

P 0»t_/_ n -i+i (j - jOKj + l)j T(/ + 1 + gig) 2 ; + 1 (r 1 - r 2 ) ; q^+ 1 fc m 

^ m _i ij ° m (21)1(21 + 1)! T(l + 2iQ) nn H ' 1 J 

We interpret i 7 ^, i 7 ^ and i 7 ^ -11 as the direct radiation from the disk, the net flux of the scattered wave caused by 
superradiant scattering, and their interference effect, respectively. To find the contribution of the scattered radiation 
in the net flux at infinity (|108|) . we pay attention to the dependence of frequency u) m . From Eqs. (|100j) . (I108[) . (|110j) 
and (1112[) . it is easy to see that in the interference effect F^ _H , the scattered radiation appear from the order of 
(Mw m ) +1 . However, this term is very smaller than the disk radiation term with the order of (Mu m )° in Eq. (|109|) . 
if the low-frequency limit is considered. 

From the viewpoint of highly energetic astrophysical phenomena, the superradiant transport of the energy flux 
to the disk will be interesting as a black hole feedback mechanism which plays a role of disk heating, while the 
contribution of the superradiant scattering of waves to the energy flux at infinity will be useful for an observational 
check of the black hole spin (see Fig. [5]). Unfortunately, as was above-mentioned, the superradiant part at infinity 
remains much smaller than the direct Mui m = Mma <C 1 is used. Therefore, it is important to analyze the case such 
that Ma ~ 1 for the frequency parameter a describing the disk radiation, by keeping the superradiance condition 
a < f^H and the regularity condition aer < 1 for the Kerr-Schild field in any region except on the disk. 

For high m modes giving m > 1 /Ma the superradiant effect may be also suppressed. Nevertheless, we must remark 
that the contribution of such modes should be taken into account if one estimates the total flux of the disk radiation 
given by the function ijj(X) and ^(X) singular at the equatorial plane. This is because as a result of the existence of 
the branch point at \X\ = 1 the coefficients a m and b m in the expansion form and ([55)1 do not rapidly decrease as 
m increase. Then, the high m modes of the vacuum field given by </>q W and </>2 W should be also efficiently generated 
from the disk radiation corresponding to the Kerr-Schild m modes. Even if the superradiant effect is small for each 
high m mode, the total sum (|61[) of the energy fluxes (£ m ) for all m modes will be an important task to discuss more 
clearly the energy transport in the disk-black hole system. Because our main purpose in this paper is focused on the 
construction of the basic formulae to calculate of the energy fluxes at the horizon, the equatorial disk and the far 
distant region, such a calculation of the total energy flux will be investigated in future works. 
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